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ABSTRACT 

We report on the development of a technique for precise radial-velocity mea- 
surements of cool stars and brown dwarfs in the near infrared. Our technique is 
analogous to the Iodine (J 2 ) absorption cell method that has proven so success- 
ful in the optical regime. We rely on telluric CH 4 absorption features to serve 
as a wavelength reference, relative to which we measure Doppler shifts of the 
CO and H 2 features in the spectra of our targets. We apply this technique 
to high-resolution (R^50,000) spectra near 2.3/im of nine L dwarfs taken with 
the Phoenix instrument on Gemini-South and demonstrate a typical precision 
of 300 m s _1 . We conduct simulations to estimate our expected precision and 
show our performance is currently limited by the signal-to-noise of our data. We 
present estimates of the rotational velocities and systemic velocities of our tar- 
gets. With our current data, we are sensitive to companions with M sin i > 2Mj 
in orbits with periods less than three days. We identify no companions in our 
current data set. Future observations with improved signal-to-noise should result 
in radial-velocity precision of 100 m s _1 for L dwarfs. 
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Introduction 



The majority of planet searches focus on sun-like stars, dwarfs of late K to early F 
spectral type. While stars smaller than these are far more common (e.g. Henry et al. 
2006), extension of the well-developed radial- velocity techniques to these cooler, dimmer 
targets poses many difficulties. Searching for planetary companions to stars with masses 
a fraction of that of the sun presents the obvious advantage that at equal orbital periods, 
companions of smaller masses can be detected with a given radial velocity precision. The 
two most successful radial-velocity techniques, the I2 absorption cell and Th-Ar comparison 
lamp methods, have both been applied to stars of early M spectral type. This has resulted in 
the d iscovery of four M- dwarf plan etary syst ems: GJ 876 (Delfosse et al.lll998at lRivera et al. 



2005b, GJ 436 (Butler et al 



20071 ). lEndl et al. 



200J), GJ 581 (IBonfils et al. 



I2OQ5J), and GJ 674 faonfils et al. 



(120031 . l2006h carried out a systematic search of 90 M dwarfs, most of 
spectral type earlier than M3. These authors achieved an average precision of 8 m s _1 and 
estimated that the rate of Jupiter-mass companions within a < 0.7AU of these objects is less 
than 1.3%, a re sult marginally inconsistent with ra dial-velocity observations of larger main 
sequence stars ( iMarcy et al.ll2005l : lUdry et al.l 120071 ). 



Ida fc Linl (120051 ) have conducted theoretical simulations of planet formation via core 



accretion, and predict that companions of Neptune-mass and smaller should be common in 
short-period orbits around lower main-sequence stars. However, their sim ulations do not 



extend dow n to the low est mass stars, or to brown dwarfs. Work by both lLaughlin et al. 



(12004] ) and iBosd (120061 ) specifically addresses the formation of planets around low-mass 
stars. They also predict that Neptune-mass companions should be relatively common, while 
Jupiter-mass companions should be much less common orbiting low-mass stars than sun- 
like stars. The microlensing detectio n of two sub-Neptune mass planets orb iting M dwarfs 
tentatively supports this hypothesis (IGould et al.l 12006c iBeaulieu et al.ll2006l ). Searching for 
giant companions to low-mass stars represents an important test of current planet formation 
theories. 

Today we know of a large number of objects populating the lowest reaches of the main 
sequence (e.g. Cruz et al. 2003a) as well as objects too low in mass to reach the main 
sequence. These objects, classified into the late M, L and T spectral types, have masses 
close to the minimum mass required for main sequence hydrogen burning. Objects in the L 
class straddle the O.O75M limit that is commonly invoked to delineate brown dwarfs from 
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stars (IBurrows et al.ll200ll ). The question of the formation of companions to these very low- 
mass stars has been little explored. There is, however, strong evidence from observations of 
infrared excesses of very low-mass stars in star-forming regions (see Luhman et al. 2007 for 
a review) that, when young, these objects possess protoplanetary disks at a rate similar to 
that of more massive stars. 

At fixed orbital period, the reflex radial velocity amplitude due to a companion of a given 
mass scales with host mass as M* 2 ^ 3 , making an L dwarf with a mass of O.1M = 105Mj a 
very accessible target for a radial-velocity search for planetary companions. Specifically, the 
semi-amplitude of the radial-velocity signal from a planet of mass M p orbiting an L dwarf 
with a period P, orbital inclination i, and eccentricity e is: 
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T he masses of fie l d L d warfs are not yet well established from purely dynamical esti- 



mates. 



Burrows et al.l (1200 f I ) find that the mass of a low-mass star may be estimated from 



its age t and effective temperature T eS according to 
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M Q ' V1550K J \ t 
If we assu me that all of our objects are older than lGyr and have T e g in the range 1800K 
to 2400K (IDahn et al.ll2002l ) we infer that all of our targets have a mass of approximately 
0.1M Q since mass depends rather weakly on T c q and age for our field (presumably old) 
objects. This may overestimate the mass of some of our targets by 15% or more, but our 
ability to detect planets only improves with lower primary mass, hence this is a conservative 
assumption. A detection of the resonance absorption line of Li in the spectrum of a low-mass 



star would constrain the mass to be less than 65 Mj (IMagazzu et al.lll993l ; iKirkpatrick et al. 



19991 ; iMartm et al.lll999l ). A literature search identified no searches for Li in any of our 
targets. Assuming M* = O.1M , a companion with M p sinz = Mj with an orbital period of 
3d would induce a radial- velocity variation with a peak-to-peak amplitude of 1.3 km s _1 . By 
optical radial velocity standards this is an enormous signal, but precision sufficient for the 
detection of this signal has never been demonstrated in observations of L dwarfs. 

The extremely low luminosity of these objects makes obtaining precise radial velocities 
difficult . iGuenther fc Wuchterll (120031 ) observed one L dwarf in the optical usin g Th-Ar com 



parison lamps as a wavelength reference and obtained a precision of 1.3 km s _1 . iBailer- Jones 
(12004] ) observed a sample of early L dwarfs in the optical and used a cross-correlation tech- 
nique to produce radial-velocity measurements at a single epoch with a precision between 
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0.5 and 4 km s _1 . The extremely red colors of these objects (V — K—10.0 at L0; Dahn et al. 
2002) motivate a move to the near infrared (NIR) and the development of a set of techniques 
that could be used to produce precise r adial velocities with spectrographs operating in that 
wavelength range. iMartm et al.l (120061 ) present observations of the young brown dwarf LP 
944-20 (spectral type M9) near 1.25/iin using Keck/NIRSPEC and report a precision be- 
tween 0.4 and 1.6 km s _1 from their cross correlation techniq ue. These authors also presen t 
optical measurements of LP 944-20, some in common with iGuenther fc Wuchterll (120031 ). 
and find that their NIR techniques are achieving similar precision to the optical techniques 
on this object. The optical radia l velocities of LP 94 4-20 exhibit significant variation while 
the NIR radial velocities do not. iMartfn et al.l (120061 ) suggest that spots, with greater con- 
trast in the optical than the NIR, are responsible for the observed radial velocity variations 
and conclude that radial-velocity measurements of low-mass stars in the optical may be 
fundamentally limited by the stars themselves. 

Here we describe a technique for high-precision radial-velocity measurements of very 
low-mass stars in the NIR and present some initial results. This technique, which uses 
the telluric absorption features of the Earth's atmosphere in place of a gas absorption cell, 
exploits the rich radial velocity content of the CO bandhead in the K-band spectra of these 
cool objects. In Section [2], we describe our observations and the extraction of the spectra. 
In Section [3l we describe our modeling procedure, and in Section 0] we present results from 
a single observing run spanning five nights, and place limits on the presence of orbiting 
companions for our targets. 



Observations 



The use of telluric absorption features as a refere nce from which to deriv e precis e radia l 
velo cities of stars has a long history in the literature. iBalthasar et al.l (119821 ). ISmithl (119821 ). 
and ICaccin et al.l (119851 ) demonstrated the use of telluric O2 with a stability of 10 m s _1 . 
The NIR region is rich with telluric features due to CH 4 . The L dwarfs contain a wide 
array of their own absorption features due to molecules such as H2O and CO. At these 
effective temperatures, the dominant carbon-bearing molecule is CO, not CH4, so there is 
no confusion between object and telluric CH4. The bandhead due to the 2 — transition 
of CO at 2.3 /im results in a set of periodic features in the L dwarf spectrum. We seek to 
exploit the superposition of the telluric CH 4 and L dwarf CO and H 2 features in order to 
extract precise radial velocities from high-resolution observations at these wavelengths. 



W e observed a sample of nine L dwarfs with the Phoenix spectrograph (IHinkle et al. 
2003J) on Gemini-South over five nights spanning UT 2006 January 10-14 under observing 
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program GS-2005B-C-3. Observing conditions were mostly free of clouds during this observ- 
ing run, and seeing remained between 0.3"and 0.8"in K band. We observed all of our targets 
at airmass between 1.0 and 1.6. We used a slit with a width of 3 pixels (0.26") to yield an 
approximate resolution of R=50,000 with the single-order spectra spanning 0.012^tm, from 
2.296/mi to 2.308/im. Exposure times ranged from 1000 to 1800 seconds per nod position 
and were scaled so that the expected S/N per pixel in the extracted spectra would be ap- 
proximately 10.0. In practice we found that the S/N of our many of our extracted spectra 
fell below this expectation by up to 50%. This is primarily due to guiding-induced jitter 
coupled to slit losses as Phoenix has no slit-viewing camera. We gathered observations in 
nodded pairs, with consecutive exposures offset in opposite directions relative to the center 
of the slit. This observing strategy facilitates subtraction of the bright sky lines. Our tar- 
gets, described in Table 1, range in spectral type from L0 and L6 and were selected from 
the sample of known L dwarfs brighter than K=12. We observed each object at least three 
times during the observing run and observed seven of the nine L dwarfs four times. We also 
observed rapidly rotating A stars each night in order to help calibrate the spectrograph PSF 
from the telluric lines superimposed upon the featureless spectra of these hot stars. 

We created a single flat field for each night by first subtracting a median dark image of 
the same exposure time from each of a nightly set of 15 individual internal flats and then using 
a median combination of normalized frames to create a final flat field image. After dividing 
each science frame by the flat field image, the sky lines and dark current were effectively 
removed from the 2-D spectra by subtracting nod pairs. We extracted the 1-D spectr a from 



the sk y-subtracted nod pairs following the optimal extraction procedure outlined by iHorne 



(119861 ) . We created a bad pixel mask through visual inspection of the flat field images in 
order to mask bad pixels during optimal extraction procedure. We propagated the flux of 
the sky lines from the individual images in order to correctly evaluate the noise estimate on 
the flux extracted from the subtracted nod pairs. We excised 200 pixels from the blue end of 
each spectrum that were contaminated by cosmetic defects and a strong sensitivity gradient 
in that region of the detector. The final 1-D spectra spanned 800 pixels centered around the 
wavelength 2.303/im. Representative spectra are plotted in Figures 1 through 3. 



3. Spectral Modeling and Fitting 



We estimated the radial velocities of the targets using a detailed modeling procedure. 
Our methods are similar to those used in the I2 cell metho d (IButler et al.lll996l ). Following 
a strategy similar to that described in IValenti et al.l (119951 ). we build up a many-parameter 
model to describe the spectrum of our target, the spectrograph response, and the Earth's 
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atmosphere. 



10 



Our model begins with two high-resolution template spectra. A high-resolution (5 x 
am pixel -1 ) spectrum of the Earth's atmosphere is provided by iLivingston &: Wallace 
(Il99ll ). This spectrum is derived from Fourier transform interferometer observations of 
the Sun at different airmasses. We also use a high-resolution synthetic spectrum of an 
L dwarf pre pared specific a lly fo r this project. These synthetic spectra are computed as 



descri bed in iMarley et al.l (120021 ). with updated opacities as described in iFreedman et al. 



( 120071 ). However, only the line opacities of CH 4 , CO, H 2 0, TiO and H 2 S are relevant in this 
wavelength region. Continuum opacities include the pressure- induced opacity of molecular 
hydrogen and helium and the Mie opacity of the iron and silicate clouds. To treat the 



clouds, the models apply the condensation cloud model of lAckerman &: Mar ley I (120011 ) with 
a sedimentation parameter of / se d = 3, corresponding to a moderate amount of condensate 
settling. Both absorption and non-isotropic scattering from the cloud s are accounte d for in 
the radiative transfer. The models used here have solar metallicity (ILoddersI 120031 ) and a 
fixed gravity of log g = 5 (cgs) and cover a range of T e g from 1200 to 2400 K. Complete 
details of the models will be presented in an upcoming publica tion. Fits of model sp ectra 
to moderate resolution data in the near-infrared are presented in lCushing et al.l (120071 ). The 
synthetic spectra provide monochromatic fluxes spaced 4.2 x 10~ 6 /mi apart. We convolve 
and re-sample the product of these high-resolution spectra to generate the model that we 
then fit to the extracted 1-D spectra. Our model has several free parameters. The parameters 
related to the L dwarf are Vsinz, T e g, and radial velocity and the parameters related to the 
spectrograph are the PSF width, flux normalization, and the wavelength solution (i.e. the 
mapping from wavelength to pixel position). 

We find that a simple Gaussian model of the PSF is sufficient for our data and addi- 
tional parameters do not significantly improve the resulting fits. We also find no evidence 
for variation of the PSF width across the spectrograph order. Such a simple model likely 
works well because the spectrograph is not cross dispersed, spans only a very small range of 
wavelength, and we have relatively few (fa 10) telluric lines, each of low S/N. We describe 
the wavelength-to-pixel mapping with a second-order polynomial function. We experimented 
with higher-order functions but found that the additional parameters did not improve the 
quality of the resulting fits. There is a strong throughput gradient in the blue-most ~ 100 
pixels of all of our spectra. Presumably due to an inadequate calibration, this was easily cor- 
rected by a two-parameter normalization process: the first parameter determines the overall 
level of the spectrum while the second parameter describes a linear, sloping normalization 
applied to the first 100 pixels. 



The first step in our model generating procedure is to Doppler shift the wavelengths 
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corresponding to each monochromatic intensity in the synthetic L dwarf spectrum. The 
synthetic L dwarf and telluric spectra are then super-sampled using cubic-spline interpolation 
onto a regular grid in wavelength at a resolution of approximately 1 x 10~ 6 /xm pixel -1 , about 
10 times the resolution of the observed spectra. We experimented with linear and quadratic 
methods of interpolation and found only small differences in the cores of absorption lines, 
which did not alter our final r esults . The L dwarf spectrum is then convolved with a rotational 



broadening kernel following iGrayl (Il992l ). We assume no limb darkening in the rotational 
broadening kernel. It is necessary to scale the telluric spectrum from airmass 2.0 to the 
airmass of each observation. We use a logarithmic scaling relation empirically determined 
from Phoenix observations of rapidly rotating A stars at a range of airmasses. The synthetic 
L dwarf and telluric spectra are multiplied together and this product is convolved with a 
Gaussian representation of the PSF of the spectrograph. The basic steps in the creation of 
the super-resolution model spectrum M(A) can be expressed 

M(A) = ( l(\x (l + ^*K xT(\)^j*PSF (3) 

where * indicates convolution, L(X) is the synthetic L dwarf spectrum, K is the rotational 
broadening kernel, T(\) is the telluric spectrum, and PSF is the spectrograph point spread 
function. 

We resample this super-resolution spectrum, which is defined in terms of wavelength, 
onto the detector grid, which is defined in terms of pixel number. This mapping is described 
by a second-order polynomial. We experimented with two methods of accomplishing this 
task. The first involved integrating the finely-sampled spectra over the extent of each detector 
pixel. This produced nearly identical results to a cubic spline interpolation for estimating 
the value at the detector position defined by the center of each pixel. We chose to use 
the cubic spline method for mapping the spectrum onto the pixel grid. Lastly, the entire 
model spectrum, now at the same resolution as the actual data, is multiplied by the two 
normalization factors described earlier. The steps in our procedure are illustrated in Figure 
1, which shows the model spectra at each stage in the modeling process. 

In order to estimate the parameters we calculate the goodness of fit between our model 
and the data and pr oduce a y 2 value . We then utilize the amoeba downhill-simplex mini- 



mization procedure (IPress et al.lll986l ) to minimize \ 2 an d estimate the best-fit parameters 
for each of our spectra. For each object we executed our fitting scheme in two phases. First, 
the V sin i value was allowed to vary and model L dwarf spectra of the 13 different effective 
temperatures were used to generate our model. We assigned to each object a T e s correspond- 
ing to the theoretical L dwarf spectrum that yielded the lowest \ 2 value. This value is not 
meant to be an estimate of the actual T e g. Using this best-fit T e g the \ 2 surface in Vsini 
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was determined for each observed spectrum by fitting the data with a range of fixed V sin i 
values while allowing the other parameters to float. A parabolic fit to this surface produced 
a Vsini estimate for each observed spectrum for each object as well as an error estimate. 
We used the weighted average of the V sini estimates from each spectrum to estimate the 
V sin i and an error set to the rms of the individual estimates. This statistical V sin i error 
likely underestimates the true error on our determination of this parameter. Experiments 
with theoretical L dwarf spectra with surface gravities other than logg = 5 indicated that 
the Vsini values we estimate depend on the surface gravity of the model. Models with 
log g = 5.5 resulted in Vsini estimates that were systematically lower by 1 — 2 km s , 
albeit with fits of significantly lower quality. This first phase of the fitting process assumes 
no degeneracy between V sini and T e g. We performed tests with artificial data designed to 
mimic the properties of our Phoenix observations (see below) and found this assumption to 
be valid. In the second phase of our fitting process we fix the values of T eff and V sin i and use 
amoeba to determine the seven remaining parameters. These seven parameters are: wave- 
length solution (3), flux normalization (2), RV(1), PSF width (1). The best-fit radial- velocity 
estimates are corrected to the solar system barycenter using the IRAF routine bcvcorr (G. 
Torres; private communication). Observed spectra and best- fit models for two of our targets 
are shown in Figures 2 and 3. 

It is not trivial to assess the accuracy with which these parameters are determined 
when using a minimization scheme like amoeba. We conducted simulations in an attempt to 
estimate the errors on our radial velocity determinations as a function of S/N and \^sin i. We 
first estimated the S/N of the extracted spectra by taking the median of the S/N values of the 
800 pixels. We then simulated Phoenix spectra of L dwarfs following the prescription above 
but with the additional step of addi ng noise, determine d by photon statistics and the gain and 



readnoise of the Phoenix detector (IHinkle et al.l 120031 ) , following the final re-sampling onto 



the pixel grid. These simulated Phoenix spectra were generated with spectrograph-related 
(PSF, wavelength solution) parameters fixed to the averages of the values determined from 
the data, at an airmass of 1.2, and radial velocity drawn from a uniform distribution between 
—20 and 20 km s _1 . We fed these simulated spectra into our amoeba parameter estimation 
code to estimate the accuracy with which the known L dwarf parameters could then be 
recovered. Here, we assumed the correct input template but experiments with different 
templates found that changing the template had little effect on the derived parameters, with 
the exception of template spectra with T c s < 1600K (below that of any of our targets) . We 
carried out the simulation over a grid of S/N from 2 to 100 and Vsini from 5 to 50 km 
s _1 . The results of this simulation are shown in Figure 4. Comparison of the scatter of the 
radial-velocity estimates of our targets to the single measurement errors predicted by our 
simulations indicates that we are achieving the expected precision. We therefore assign errors 
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to our radial- velocity measurements based on interpolations of the results of our simulations. 



Results 



For each of our targets we find rotation velo c ities i n exc ess of the 6 km s 1 velo city 
resolution of our spectrograph. Like iBailer-Joned (120041 ) and iMohantv fc Basri J2003h we 



find th at our L dwarfs are relatively rapid rotators compared to field M dwarfs (IDelfosse et al. 



1998bl ). We estimate Vsini values between 7.06 and 20.68 km s" 1 for our nine targets. We 
present the best fit Vsini values and systemic velocities (V 7 ) for each of our targets in Table 
1 . The errors on V sin i are derived from the fits described in Section [3] and are based on 
the standard deviation of the V sin % values from the individual observations. The value for 
Vy is the average of the individual radial-velocity measurements, and the error we quote is 
the average error divided by N 1 ^ 2 , where N is the number of individual observations. The 
individual measured radial velocities are presented in Table 2 and the radial velocity curves 
for each target are shown in Figures 5,6, and 7. For six of our targets we achieve a radial 
velocity precision less than 450 m s -1 per individual measurement. With our typical S/N 
between 5 and 15, it is clear from Figure 4 that we are currently limited by the noise of 
our observations and not by systematic effects. Increasing our S/N to 20 would improve our 
radial velocity precision to 150m s^ 1 . As a result, avenues for future improvement in our 
precision are promising. 

None of the objects show evidence for radial velocity variability. In each case the fit 
to the data of the null hypothesis of constant radial velocity is good with xl ~ 1-0- There 
is evidence that the error estimates derived from our Monte-Carlo simulations are modestly 
overestimated, as three of our targets have xl < 0.5. The source of this discrepancy is not 
immediately clear. An incorrect estimate of either Vsini or the median S/N of the 800 
pixels in the spectra could lead to an incorrect estimate of the errors. We note, however, 
that these low x 2 values are not unexpected considering the small number of observations 
of each object and the total number of objects observed. Our lowest x 2 value (2M1048+01, 
X 2 = 0.34) has a probability of 6.4%, which is not surprising. 

While we don't have a sufficient number of targets in our current sample to place 
interesting limits on the frequency of close-in, Jupiter-mass companions to L dwarfs, we can 
quantify the masses of companions that we would have securely detected in the observations 
of our nine targets. We carried out such a simulation by injecting sinusoidal radial-velocity 
signals into the data for each target and examining the statistics of the resulting radial- 
velocity curves. For each target we created 10 3 fake data sets by the addition of a signal 
due to a planetary companion with 0.5 < Msinz < 10Mj, a period 1.0 < P < 3.0 days, 
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and random orbital phase. For each simulation we assume an L dwarf mass of 100 Mj 
and estimate the value of Msinz at which we would be 95% likely to detect a deviation 
in x 2 with a confidence of 95%. We find that for all of our targets we can effectively 
rule out companions of Msinz > 4.5Mj with periods of less than 3 days and that for 
four of our targets (2M0109-51, 2M0835-08, 2M0355+11, 2M1045-01) companions with 
Msinz > 2.0Mj can be excluded in that period range. 



5. Conclusions 



We describe a method to estimate the multiepoch radial velocities of L dwarfs with 
a precision sufficient to detect short-period Jupiter-mass companions. Our technique uses 
atmospheric telluric CH4 absorption features as the wavelength reference, and the molecular 
CO and H 2 absorption spectrum in the L dwarf to estimate the radial-velocity offset. We 
have demonstrated an average precision of 300 m s" 1 over a time span of five nights. Monte- 
Carlo simulations show that even a modest increases in the S/N of our data should lead to 
a factor of 3 improvement in our precision. This level of precision opens up the possibility 
of detecting Jupiter-mass planets with orbital periods of months. Improvements in S/N and 
radial-velocity precision can be realized with current and future higher-throughput instru- 
ments on large telescopes. An understanding of the frequency of planetary companions to 
objects at the bottom of the main sequence will provide important insight into the process of 
planet formation. Moreover, the indirect detection of such planets orbiting nearby low-mass 
stars and brown dwarfs affords the opportunity for direct study of the planetary emission 
and characterization of the atmosphere. The ratio of the infrared brightness of a giant planet 
to its host star (smaller than 10~ 5 at A < 10/im for a giant planet orbiting sun-like star at 
0.5AU; Burrows et al. 2004) is improved by several orders of magnitude when the host star 
is an L dwarf. If it is possible to overcome the overa ll faintness of the system, direct de- 
tection of the planet through high- resolution imaging (IBeuzit et al.l 120071 ) may be possible. 
Furthermore, if the planet is found to transit its h ost, observations of secondary eclipses 
( IDeming et al.l 120051 . 120061 ; ICharbonneau et al.l 120051 ) will result in direct measurements of 
the planet and its properties. 
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Fig. 1. — Illustrative example of the steps in our modeling process: (a) is the rotationally 
broadened L dwarf theoretical spectrum , (b) is the observed KPNO telluric spectrum de- 
scribed by iLivingston fc Wallace! (1199 ll ). (c) is the product of (a) and (b), and (d) is the 
spectrum (c) convolved with the PSF. The data (small points) in spectrum (d) are the 
observed spectrum from our Phoenix observations of 2M1155-37. 
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Fig. 2. — Example of a model fit to an observed spectrum of 2M0109-51. The top panel 
shows the observed data as small points with the best-fit model (determined by the amoeba 
fitting procedure) overplotted as a line. The bottom panel shows the residuals of this fit. 
The error bar in the upper left corner approximates the expected noise (photon and read 
noise) of the spectra. We estimate the V sinz of 2M0109-51 to be 14.58±0.52 km s _1 . 
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Fig. 3. — Example of a model fit to an observed spectrum of 2M 1045-01. The top panel 
shows the observed data as small points with the best-fit model (determined by the amoeba 
fitting procedure) overplotted as a line. The bottom panel shows the residuals of this fit. 
The error bar in the upper left corner approximates the expected noise (photon and read 
noise) of the spectra. We estimate the Vsini of 2M1045-01 to be 7.06±0.71 km s _1 . 
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S/N 

Fig. 4. — Results from our Monte-Carlo simulation to estimate the radial-velocity precision 
as a function of median S/N per pixel and Tasini Each line represents a different Vsini 
value in km s _1 . Our data have S/N between 5 and 15 and Vsini between 7 and 21 km 
s -1 , corresponding to the shaded region of the plot. The RMS of the radial velocities of 
each target, along with the average S/N of the observations, are plotted as solid squares. In 
general, our achieved precision follows the expectations from our simulations. The RMS of 
the observations of 2M0109— 51 is about 150 m s -1 larger than expected. This large RMS 
results from the pair of observations taken at the third epoch (see Figure 5). 
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Fig. 5. — Radial- velocity measurements from each position of nodded pairs of spectra derived 
from the modeling processes described in Section [3l The systemic velocity V 1 has been 
subtracted. A barycentric correction has been applied and the errors are estimated from the 
results of the Monte-Carlo simulations summarized in Figure 4. 
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Fig. 6. — Radial- velocity measurements from each position of nodded pairs of spectra derived 
from the modeling processes described in Section [3l The systemic velocity V 1 has been 
subtracted. A barycentric correction has been applied and the errors are estimated from the 
results of the Monte-Carlo simulations summarized in Figure 4. 
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Fig. 7. — Radial velocity measurements from each position of nodded pairs of spectra derived 
from the modeling processes described in Section [3l The systemic velocity V 1 has been 
subtracted. A barycentric correction has been applied and the errors are estimated from the 
results of the Monte-Carlo simulations summarized in Figure 4. 
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Table 1. Sample 



Object 


RA 


Dec 




Sp. Typ. 


Template T a g 


V sini 


V 7 




rrns(RV) 


S/N 




J2000 


J2000 






K 


km s 


km s 1 




m s 




2M0109-51 


01:09:02 


-51:00:49 


11.10 


L2 a 


1900 


14.58+0.52 


-1.30+0.20 


1.51 


570 


9.6 


2M0306-36 


03:06:12 


-36:47:53 


10.63 


L0 a 


1900 


20.68+0.65 


+11.58+0.09 


1.15 


260 


14.2 


2M0314+16 


03:14:03 


+16:03:06 


11.24 


L0 b 


2000 


17.72+1.32 


-6.80+0.24 


0.74 


640 


5.5 


2M0355+11 


03:55:23 


+11:33:44 


11.53 


L6 b 


1900 


10.31+0.53 


+12.24+0.18 


0.89 


440 


7.9 


2M0523-14 


05:23:38 


-14:03:02 


11.64 


L5 C 


2100 


15.26+0.60 


+11.82+0.16 


0.63 


400 


8.4 


2M0835-08 


08:35:43 


-08:19:24 


11.14 


L5 d 


1600 


16.92+0.50 


+29.96+0.14 


0.63 


390 


9.8 


2M1045-01 


10:45:24 


-01:49:58 


11.78 


Ll c 


2200 


7.06+0.71 


+6.33+0.18 


1.56 


510 


7.0 


2M1048+01 


10:48:43 


+01:11:58 


11.62 


L4 f 


2200 


10.29+0.47 


+24.09+0.10 


0.34 


280 


8.2 


2M1155-37 


11:55:40 


-37:27:35 


11.46 


L2 C 


1600 


15.93+0.67 


+45.63+0.13 


0.50 


360 


8.7 



Note. — The first five columns are object-specific information gathered from dwarfarchives.org. The reported spectral types are based 
on NIR observations when available, optical observations otherwise. The values of T e fj represent not an estimate of the actual effective 
temperature of the object but the temperature of the model spectrum that best fit the data and that was used in the radial velocity 
analysis. Values of x 2 are per degree of freedom and represent the fit to the null hypothesis of constant radial velocity. The values of 
rms(RV) represent the standard deviation of the radial-velocity measurements for each object treating each nodded position as a separate 
observation. The values of S/N are the average of the median signal-to-noise per pixel of the individual observations of each object. The 
errors on Vsini and V 7 are statistical and do not take into account likely systematic errors. 

HLodieu et ail l l2005f) 

ISchmidt etafl ll2007ll 

IWilsonl J2003h 

fcruz et al.1 d2003bh 

jGizisI j2002h 

^Kendall et al] lb004T) 
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Table 2. Radial Velocity Measurements 



Object 


HJD-2400000 


RV 


a 






km s~ 1 


km s _1 


2M0109-51 


53744.54448 


-1.25 


0.31 


2M0109-51 


53744.55882 


-1.39 


0.34 


2M0109-51 


53745.53691 


-1.43 


0.43 


2M0109-51 


53745.55125 


-1.79 


0.46 


2M0109-51 


53747.52939 


-2.33 


0.39 


2M0109-51 


53747.54720 


-0.69 


0.38 


2M0109-51 


53748.53308 


-0.80 


0.46 


2M0109-51 


53748.54742 


-0.75 


0.50 


2M0306-36 


53745.65385 


12.01 


0.28 


2M0306-36 


53745.66819 


11.42 


0.28 


2M0306-36 


53746.62129 


11.70 


0.21 


2M0306-36 


53746.63330 


11.23 


0.19 


2M0306-36 


53747.62901 


11.29 


0.45 


2M0306-36 


53747.64104 


11.64 


0.42 


2M0306-36 


53748.61104 


11.59 


0.23 


2M0306-36 


53748.62537 


11.78 


0.27 


2M0314+16 


53744.58707 


-7.47 


0.70 


2M0314+16 


53745.58148 


-5.94 


0.63 


2M0314+16 


53745.59583 


-6.23 


0.57 


2M0314+16 


53746.58386 


-6.79 


0.59 


2M0314+16 


53746.59820 


-7.69 


0.68 


2M0314+16 


53748.57643 


-6.93 


0.66 


2M0314+16 


53748.59077 


-7.23 


0.64 


2M0355+11 


53745.61973 


12.36 


0.72 


2M0355+11 


53745.63407 


12.82 


0.83 


2M0355+11 


53746.53400 


12.11 


0.43 


2M0355+11 


53746.55528 


12.48 


0.54 


2M0355+11 


53747.58124 


12.16 


0.48 


2M0355+11 


53747.60254 


11.51 


0.49 


2M0523-14 


53746.65617 


12.31 


0.37 


2M0523-14 


53746.67747 


11.53 


0.42 


2M0523-14 


53747.66353 


12.02 


0.47 


2M0523-14 


53747.68134 


11.69 


0.48 


2M0523-14 


53748.65352 


12.10 


0.41 


2M0523-14 


53748.67480 


11.25 


0.56 


2M0835-08 


53744.68678 


30.33 


0.56 


2M0835-08 


53744.70112 


30.05 


0.55 


2M0835-08 


53746.71119 


29.94 


0.29 


2M0835-08 


53746.72554 


30.15 


0.30 


2M0835-08 


53747.70936 


29.21 


0.46 


2M0835-08 


53747.72371 


29.64 


0.49 


2M0835-08 


53748.70701 


29.96 


0.48 


2M0835-08 


53748.72484 


30.42 


0.62 


2M1045-01 


53744.74669 


6.04 


0.64 


2M1045-01 


53744.76797 


7.07 


0.59 
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Table 2 — Continued 



Object 


HJD-2400000 


RV 


a 








km. s 


2M1045-01 


53745.77659 


6.23 


0.82 


2M1045-01 


53745.79789 


5.60 


0.81 


2M1045-01 


53746.76673 


6.01 


0.44 


2M1045-01 


53746.78455 


6.92 


0.54 


2M1045-01 


53747.80330 


6.70 


0.57 


2M1045-01 


53747.81996 


6.10 


0.72 


2M1048+01 


53744.79635 


23.85 


0.41 


2M1048+01 


53744.81070 


23.55 


0.43 


2M1048+01 


53746.81080 


23.98 


0.36 


2M1048+01 


53746.82747 


24.19 


0.37 


2M1048+01 


53747.76434 


24.32 


0.37 


2M1048+01 


53747.77984 


24.13 


0.40 


2M1048+01 


53748.77484 


24.31 


0.40 


2M1048+01 


53748.79613 


24.36 


0.37 


2M1155-37 


53744.83327 


45.13 


0.46 


2M1155-37 


53744.84647 


45.55 


0.44 


2M1155-37 


53745.82581 


45.24 


0.43 


2M1155-37 


53745.84248 


45.38 


0.46 


2M1155-37 


53746.85162 


46.01 


0.40 


2M1155-37 


53746.86773 


45.84 


0.47 


2M1155-37 


53747.84184 


45.86 


0.46 


2M1155-37 


53747.85617 


46.06 


0.44 



Note. — Individual radial-velocity observations. 
Barycentric corrections have been applied to the listed 
here. The errors on each measurement are estimated 
from the Monte-Carlo simulations described in Section 
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